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Abstract 

The large reservoir of antibiotic resistant bacteria in raw and treated water supplies is a 
matter of public health concern. Currently, the National Antimicrobial Resistance Monitor- 
ing Systems, a collaborative effort of the Centers for Disease Control, the US Department 
of Agriculture, and the US Food and Drug Administration, does not monitor antimicrobial 
resistance in surface waters. Given the serious nature of antibiotic resistance in clinical 
settings, and the likelihood that antibiotic resistant bacteria can be transmitted to humans 
from large environmental reservoirs via drinking water, explanations for the distribution of 
antibiotic resistant bacteria and tools for studying this distribution must be found. Here we 
focus on mathematical modeling of cultivable bacteria in a river, which will be used to study 
the distribution of antibiotic resistant bacteria in the environment. We consider both an- 
tibiotic resistant and non-antibiotic resistant bacteria in the model, and, taking into account 
the strong correlation between land use and antibiotic resistant bacteria in rivers, we include 
a function for the influx of bacteria into the river from the shore. We simulate the model for 
two different time scales and show that if there is too many bacteria from the land entering 
the river, the river entirely fills with antibiotic resistant bacteria, while less frequent influxes 
allows time for the bacteria to lose the antibiotic resistant gene. This mathematically verifies 
that reduction in antibiotic use near the banks of rivers, will reduce the counts of antibiotic 
resistant bacteria in rivers. 

Keywords: mathematical model, non-linear population dynamics, antibiotic resistance, 
bacteria, riverine system 



1. Introduction 

The age of antibiotics is usually traced back to 1928 - the year that Alexander Fleming 



discovered penicillin. The first clinical uses of penicillin occurred in the early 1940's (Dawson 



et al. , 1941 ; Parascandola, 1997; Grossman , 2008 ) , and penicillin resistance was first reported 
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in hospital isolates in 1947, only a few years after its introduction (Lewis, 1995). Since 
that time, antibiotic resistance in clinically relevant bacteria has become a major health 
concern. The US Food and Drug Administration estimates that approximately 70 percent 
of pathogenic bacteria encountered in hospital settings are resistant to one or more of the 



drugs commonly used in their treatment (Bren, 2002). Concern over the spread of antibiotic 



resistant bacteria led to the formation of the National Antimicrobial Resistance Monitoring 



System (NARMS) in 1997 ( [Bartholomew et al.[ |2006[ ). NARMS is a collaborative effort of 
the Centers for Disease Control, the US Department of Agriculture, and the US Food and 
Drug Administration. Although NARMS monitors antibiotic resistant bacteria isolated from 
humans, food animals, and raw food products, it does not monitor antimicrobial resistance 
in surface waters. 

The nation's surface waters are vital, multiple-use resources. Major rivers, in particular, 
are used as transportation routes, recreational venues, and industrial water sources, as well 
as, receptacles of stormwater runoff, sanitary sewage, and industrial wastes. In addition, sur- 
face waters are the sources of drinking water for millions of people worldwide. For example, 
the Ohio River Valley Water Sanitation Commission estimates that the Ohio River provides 



drinking water for more than five million people (ORSANCO, 2008). Several studies have 



shown that large numbers of freshwater bacteria found in these life sustaining waterways 



are resistant to one or more antibiotics (Cooke, 1976 Ash et al. , 2002 Niemi et al. 1983 



Kelch and Lee, 1978 Bennett and Kramer, 1999 Smith, 2006; Dotson, 2008). The resis 



tance found in freshwater bacteria is problematic for humans since, it has been shown that 



[ antibiotic resistant bacteria can be isolated from tap water (Armstrong et al. , 1981 Nagy 



et al. , 2009). Therefore, the large reservoir of antibiotic resistant bacteria in raw and treated 



water supplies is a matter of public health concern. 

Many studies have correlated the presence of antibiotic resistant bacteria in surface waters 
with the human activities occurring along the shore. Land uses that have been correlated 
with antibiotic resistance include wastewater treatment plants (Frank, 2005), urban areas 



and heavy metal pollution (McArthur and Tuckfield, 2000), and agricultural pastures and 



row crops (Kuhn et al. , 2005 



and other impervious surfaces (Boon and Cattanach, 1999 Somerville et al., 2004), industrial 



Dotson, 2008). As an example, regarding the Mud River, WV, 



Dotson specifically studied the watershed land use in relation to the count and distribution 
of antibiotic resistant bacteria, and concludes "the evidence shows a correlation of antibiotic 



resistant bacteria to areas of livestock in the watershed (Dotson, 2008)." 



The source of environmental antibiotic resistant strains has been attributed to the use 
of antibiotics both in medical and agricultural settings. Therefore, previous studies have 
assumed that aquatic antibiotic resistance derives directly from human and animal fecal 



contamination entering surface waters via inadequate sewage treatment or runoff (Cooke 



1976 Ash et al. , 2002 Raloff, 1999). If this assumption holds true, then antibiotic resis- 
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tant bacteria in an environmental sample should be a subset of fecal bacteria, and their 
distribution should be predicted using standard methods for the microbiological assessment 
of water and wastewater. However, recent studies have shown that cultivable antibiotic 
resistant bacteria are not a subset of fecal indicator bacteria, which implies that their dis- 



tribution in freshwater environments cannot be predicted by standard techniques (Smith 



2006 Somerville et al. 2007; Dotson, 2008) 



We conclude, there is a large population of antibiotic resistant bacteria in surface waters 
and drinking water that are not monitored under NARMS surveillance, and whose distribu- 
tion cannot be predicted by standard water quality analyses. Given the serious nature of 
antibiotic resistance in clinical settings, and the likelihood that antibiotic resistant bacteria 
can be transmitted to humans from large environmental reservoirs via drinking water, tools 
for modeling the distribution of antibiotic resistant bacteria in the environment are needed. 



Mahloch (1974) offers one of the first comparisons of mathematical models of bacteria 



in rivers. Mahloch compares six models for the count and distribution of coliform bacteria 
in rivers, three deterministic and three stochastic, with data collected from the Leaf River, 
MS. The six models do consider influx of bacteria from tributaries, though not from other 
sources, such as runoff from agricultural pastures and row crops. Also, the models are for 
counts of bacteria only, not for bacteria with or without an antibiotic resistant gene. 

Since Mahloch, models have been introduced to study the transfer of the antibiotic re- 
sistant gene, both at the bacteria level and at the human level. The mechanisms of gene 
transfer have been modeled with a mass action term (Levin et al. 1979) and later with 



a Michaelis-Menten kinetic (Andrup et al. , 1998). Transfer of antibiotic resistance among 



humans is most notable in hospitals. In the hospital setting, Webb et al. (2005) present 



a model which captures both the patient level and the bacteria level, in the sense of the 
bacterial load contributed by each patient. We know of no models, of the sort we present 
here, that focus on antibiotic resistant bacteria in a river with the possibility of bacteria 
entering the river due to human activity along the shore line. 

In this paper we present a novel model for cultivable bacteria in a river focusing on 
the influx of bacteria due to human activities along the shore. We consider both antibiotic 
resistant and non-antibiotic resistant bacteria in the model. Taking into account the strong 
correlation between land use and antibiotic resistant bacteria in water, we include a function 
for the influx of bacteria into the river from the shore. Our model includes the possibility 
of both vertical and horizontal transfer of the antibiotic resistant gene, and the fact that 
the land bacteria are not adapted to the river ecosystem and will not survive in the river. 
In Section [2] we describe the model and the model implementation. Realistic parameters 
are discussed and the model is simulated in Section El We conclude with our discussion in 
Section HI 



2. Model assumptions and implementation 

Throughout our model description we use the phrase "antibiotic resistant" to mean re- 
sistance in the clinical sense to a particular unspecified antibiotic, such as, tetracycline, and 
"non-antibiotic resistant" to mean not resistant to that particular antibiotic. In general, any 
one bacterium could be non-resistant to all antibiotics, or resistant to any number of antibi- 
otics at once. Though we consider only clinical resistance to one particular antibiotic, this 
model could easily be extended to consider bacteria with multiple resistance or to compare 
clinical and weaker forms of resistance. 

We use the independent variable time t in our model to represent time from the head of 
the river. As time passes, we assume the bacteria are transported down the river in such 
a way that all bacteria in the same cohort at time t stay together as they travel down the 



99 river. Typical data collected from rivers (e.g. Dotson (2008); Smith (2006)) mark the water 

ioo sample collection points using mile marks along the direction of the river, not time, which 

ioi leads one to believe that distance should be the independent variable. However, typically 

102 mathematical models use time as the independent variable, as we do here. This should cause 

103 the reader no concern, since by using river flow rates, distance from the head of the river 

104 can be converted to time from the head of the river. 

105 Only some bacteria are able to be cultivated and counted using current methods, and 
we so, throughout, we consider only cultivable bacteria in the river. In general, there are 
107 many more bacteria in a river that cannot be cultivated and counted. The main model 
loa assumptions are schematically represented in Figure [TJ Details of the model assumptions 
109 and their mathematical implementation are presented below. 

no We consider two distinct classes of cultivable bacteria in the river: those that are "river" 

in bacteria, R, and those that are "land" bacteria, L. The river bacteria are naturally occurring 

n2 bacteria which are adapted to survive in the river. The land bacteria enter the river from the 

n3 shore and we assume they are not adaptable to the river. Mathematically, the land bacteria 

lw have an additional death rate, which is large enough to offset any reproduction. Therefore, 

us the only increase in the land bacteria in the river comes in the form of immigration from the 

lie shore. In the river, all bacteria are homogeneously mixed, so any bacterium can come into 

H7 contact with any other bacterium. We assume that neither type of bacteria can transform 

us into the other. 
n9 Both the river and land bacteria are further subdivided into those that have the antibiotic 

120 resistant gene, called "resistant", Ri, Lj, and those that do not, called "susceptible", Rs, L$. 

121 We always have one particular gene in mine, for example, the gene to make a bacterium 

122 resistant to tetracycline. In this model, we are most interested in assessing human activities 

123 along the shore in relation to the antibiotic resistance in rivers, so we assume the land 

124 bacteria which are antibiotic resistant gained the gene through human activity, though some 



125 may naturally be resistant (see for example D'Costa et al. (2006)). 



126 The bacteria with the antibiotic resistant gene can transfer the gene to the susceptible 

127 bacteria and the bacteria with the antibiotic resistant gene can lose the gene. A resistant 

128 bacterium of either type can transfer the gene to a susceptible bacterium of either type. 

129 The transfer of the antibiotic resistant gene is similar to the transfer of a disease: one 

130 'infected' (resistant) bacterium comes into contact with a 'susceptible' bacterium and with 

131 some probability the result of the contact is two 'infected' (resistant) bacteria. In particular, 

132 a resistant bacterium does not lose the antibiotic resistant gene when the gene is transferred. 

133 An susceptible-infected (SI) disease model is used to represent the transfer of the antibiotic 



134 resistant gene. Gonzalo et al. (1989) shows that antibiotic resistant bacteria survive less 

135 in less polluted water. Therefore, we assume that the loss of the antibiotic resistant gene 

136 depends on the pollution, P, in the river. The loss rate is B(P), where the function has 

137 the property that as P increases B decreases and as P decreases B increases. Many such 

138 functions exists, and we use in this model B(P) = , so that pollution is measured 

v ' Ls + Lj + V 

139 by the bacteria in the river from the shore. 

wo Both the river bacteria and the land bacteria compete for space resources in the river 

141 leading to a river carrying capacity. The carrying capacity for the river used in the model 

142 is the carrying capacity for only the cultivable bacteria; in general, a river will have a larger 

143 carrying capacity for all bacteria including those that are not cultivable. We assume that 

144 the river has a constant carrying capacity for the length relevant to these interactions. A 

145 logistic growth model is used to represent the populations' dependence on the river carrying 
we capacity. The river carrying capacity can allow an increase in the number of river bacteria, 

147 but not in the number of land bacteria. We assume that the additional death rate of the 

148 land bacteria is large enough to offset any population growth that occurs due to the river 

149 carrying capacity. 

150 Since we assume that neither river nor land bacteria can transform into the other, the 

151 influence of the land bacteria (river bacteria) on the river bacteria (land bacteria) occurs in 

152 two ways: through contact and transfer of the antibiotic resistant gene, and through total 

153 population size in accordance with the river carrying capacity. 

154 Land bacteria can enter the river at time t. The increase in land bacteria in the river 

155 is split into those with the antibiotic resistant gene and those without, Fj(t) and Fg(t), 

156 respectively. We assume that Fj(t) > and Fs(t) > for all t. Generally, these rates are 

157 not constant and depend on time (location down river). Mathematically, these terms provide 

158 an external forcing term for the model. 

159 For simplicity, we assume that all rates which effect both river and land bacteria are 
wo constant and are the same for both types of bacteria. For example, the antibiotic resistance 
lei gene transfer rate is a constant a for both types of bacteria. 



162 Let the total cultivable bacteria in the river be denoted TCB = R$ + Ri + L$ + Lj. 

163 Note that TCB is not a constant since bacteria are added to the river from the land. Set 

164 the pollution measure to be L = L$ + Tj. Our model can be described with the following 

165 system of differential equations. 

<£ = - Q * s(M + (^) fl , +r (i-^) fls 

^ = F s (t)- 1 L s -aL s (R I + L I )+U^L I +r(l-'^p\L s 
§ = m - lLl + aLs(Rl + Ll) -^y i + r(l- T -^)L, 

166 The parameter interpretations are: a is the transfer rate of the antibiotic resistant gene; /3 

167 is the loss rate of the antibiotic resistant gene; r is the demographic rate due to the river 

168 carrying capacity; Fs(t) and Fj(t) are the rates of bacteria entering the river from the shore; 

169 and 7 is the death rate due to the land bacteria not being adapted to survive in the river. 
no The parameter interpretations, along with their dimensions, are also presented in Table [TJ 

171 (In Table III AR is used as an abbreviation for antibiotic resistant). We assume that the land 

172 bacteria death rate 7 is large enough to offset any population growth that occurs due to the 

173 river carrying capacity; specifically, when 7 > r any reproduction, mathematically due to 

174 the carrying capacity term, will be offset by the bacterial death due to not being adapted to 

175 survive in the river. 

176 It is easy to verify that the non-negative state space is invariant, that is, any initial 

177 condition (a set of four population levels) that is non-negative results in population levels 

178 that are non-negative for all time. 

179 3. Model Simulation 

wo The mathematical model presented here can qualitatively and quantitatively describe the 

i8i distribution of antibiotic resistant bacteria in a river effected by bacteria entering the river 

182 from the shore, and the subsequent transfer and loss of the antibiotic resistant gene. 

183 We demonstrate the model with a simulation of antibiotic resistant and non-antibiotic re- 

184 sistant land bacteria entering a river which initially has no antibiotic resistant bacteria. Most 

185 model parameters are determined from literature on survival and reproduction of bacteria in 

186 rivers, and on transfer and loss of antibiotic resistant genes in rivers. The parameter values 
used in our simulations and their corresponding references are given Table [2} For parameters 



187 



6 



r and 7 the values were determined using information on the half-life of bacteria in rivers 
found in Hendricks (1972) and the model developed for the Ohio River by LimnoTech| (2001) 



respectively. The transfer rate a applies in situations when the water (river specifically) is 



dilute; it was taken from Andrup and Anderson (1999) who assume that time is measured 



in hours. The antibiotic resistance loss rate is known to be higher in less polluted water 



193 (Gonzalo et al. , 1989). Thus we set the loss rate to be a function of pollution, in the model 

194 pollution is represented by the land bacteria. If the pollution reaches 0, then the loss rate 

195 is = 100, while if the pollution is at the river carrying capacity, then the loss rate is very 
we small, 0.0001. We assume a river carrying capacity of 1,000,000. Before the influx of any 

197 land bacteria, the river is assume to contain only non-antibiotic resistant bacteria and these 

198 bacteria have reached the carrying capacity, that is, -Rs(O) = K, Ri(0) = 0, Ls(0) = 0, and 

199 L/(0) = 0. 

200 To simulate the influx of land bacteria into the river, we assume that the bacteria enter 

201 the river from the shore at four distinct times, evenly spaced over the length of the river. 

202 At each of these influx times, the rate of influx of antibiotic resistant bacteria was chosen 

203 randomly to be between and 25. The rate of influx of non-antibiotic resistant bacteria was 

204 also chosen randomly to be between and 50, with the condition that the rate of influx of 

205 non-antibiotic resistant bacteria is larger than the influx of antibiotic resistant bacteria. The 



fact that Fs > Fi corresponds with data collected from the Mud River, WV ( |Dotson , 2008). 



At all other integer points, the influx is set to zero. The functions Fg(t) and Fj(t) are linear 
interpolations of these point sources. 

Two different simulations are presented in Figures [2] and [3j In both of these simulations, 
the same bacteria influx functions are used, while two different time intervals are considered. 
The first has total time 100 and the second has time 1,000. In the first case, there is too 
much pollution and, if the simulation were to continue, all of the bacteria in the river will 
eventually have the antibiotic resistance gene. Most of the bacteria will be river bacteria, 
Ri, though some will be land bacteria, Lj. In the second case, the river has enough time 
between the influx times for most antibiotic resistant bacteria to lose the antibiotic resistant 
gene. Between each influx, the river recovers and most bacteria are non-antibiotic resistant 
river bacteria, Rs- 

Using these parameters and influx functions, and both of the time scales, a local parame- 
ter sensitivity analysis was performed for the model. Each parameter was adjusted by ±20% 
(the functions Fs and Fj were considered together as one parameter and were adjusted up 
or down together) and the effect of this change was studied graphically. For all populations 
and both time scales, the river carrying capacity rate, r, is the least influential parameter 
of the model, having no effect on the population levels throughout the entire simulation(s). 
Each of the other parameters, a, /3, K, 7, and the functions Fs and Fj, was more influential 
on the longer time scale 1000 than on the shorter time scale 100. 

7 



For the shorter time scale, the number of non-antibiotic resistant river bacteria, R$, 



j, mv j.j.u.j.j.iu^i vi "vii uimiuivuiv i^uikmuiiu aiv^a uuvuviiu, -"'iS) 



227 was overall the least influenced population, the population tends to in every case, while 

228 the number of non-antibiotic resistant land bacteria was the most influenced population. 

229 Figure [4] shows the effect of the parameter sensitivity analysis on the non-antibiotic resistant 

230 land bacteria, L$, on the shorter time scale. For the longer time scale, R$, Ri, and L$ 

231 were very influenced by changes in all parameters, except r, while Lj was only influenced by 

232 changes in 7 and the influx functions. Figure [5] shows the effect of the parameter sensitivity 

233 analysis on the non-antibiotic resistant river bacteria, Rs, on the longer time scale. 

234 4. Discussion 

235 The model presented here is, by its very definition, a generalization of the actual processes 

236 that effect the amount and distribution of antibiotic resistant cultivable bacteria in rivers. We 

237 focus here on three aspects of antibiotic resistant bacteria in rivers: the influx of bacteria into 

238 the river from the shore, the transfer and loss of the antibiotic resistant gene, and the river 

239 carrying capacity. The model presented captures these aspects well and, by making some 

240 simplifying assumptions, the model uses only five parameters and the unknown functions F$ 

241 and Fj. 

242 The model simulations demonstrate that human involvement can have a significant influ- 

243 ence on the level of antibiotic resistance of bacteria in rivers. After each of the land bacteria 

244 influx points, the level of antibiotic resistant river bacteria increases, while the non-antibiotic 

245 resistant river bacteria decreases. And, only with sufficient time between influxes, can the 

246 bacteria lose the antibiotic resistant gene and the river recover. Sources of antibiotic resis- 

247 tance entering rivers must be determined and the corresponding activity adjusted to prevent 

248 the resistance gene from becoming typical in a river. Without such interventions, antibi- 

249 otic resistance will spread and clinical uses of antibiotics may become ineffective, leading to 

250 serious public heath crises. 

251 There are aspects of the transfer of the antibiotic resistant gene between bacteria that 

252 are not captured by this model. For example, we assume (implicitly) that when a bacterium 

253 'dies' it can no longer transmit the antibiotic resistant gene to another bacterium. In reality, 

254 some bacterial death can cause pieces of a bacterium's cell makeup to disperse into the river, 

255 where it is possible for some other organisms to collect these dispersed pieces, including 

256 the antibiotic resistant gene. For another example, some researchers distinguish between 

257 antibiotic resistant (in the clinical sense) and weakly antibiotic resistant; this distinction 

258 does not appear in the model presented here, though extra population classes would be easy 

259 to add to the model for to separate clinical resistance from weakly resistant. Despite its 

260 limitations, we believe that the model captures some crucial aspects regarding the dynamics 



of antibiotic resistant bacteria in rivers, specifically involving human use of antibiotics near 
rivers. 

In this paper, we offer a model which is only a preliminary formalism of the river and 
antibiotic resistance situation, however, we believe that the main aspects of the situation have 
been accounted for in this model. This novel model opens a wide range of possible future 
directions for ecological study. The model and model simulations suggest some insights 
and refinements of the model, which will lead to better understanding of the count and 
distribution of antibiotic resistant bacteria in rivers. 

1. Due to the uncertainty of the influx of the land bacteria and the sensitivity of the 
model to the influx functions Fs(t) and Fj(t), care must be taken when these functions 
are determined. Studies have shown that activity along the shore, such as farming 
row crops, livestock pastures, and the presence of urban areas, has a strong positive 



correlation with the amount of antibiotic resistant bacteria are in the river (Dotson 



2008). Assessing this dependence and using the relationship to generate accurate influx 
functions needs to be addressed. 

2. Most rivers changes size, both width across and depth, throughout their length, which 
leads to the conclusion that the carrying capacity should be a function of the length 
of the river, time t in this model. For a given river, some measure of carrying capacity 
must be found and measured along the length of the river, and used to create a river 
carrying capacity function K(t) that depends on location along the river. 

3. One obvious refinement would be to include spatial aspects in the model. Here we 
assume, basically, that the river is one dimensional; it just has a length, which cor- 
responds to time in the model. In the river, bacteria homogeneously mix, and all 
bacteria in the cohort at time t travel down the river together. From the shore, the 
land bacteria immediately enter the river and mix with the bacteria already present. 
A more realistic model would include the spatial spread of the bacteria from the shore, 
and the possibility of spreading down the river at different rates, not traveling together 
as a cohort. 

We assert that this model could be generalized to any situation with non-native organisms 
entering an environment that they are not adapted to live in and their subsequent interactions 
with the native organisms. The transfer and loss mechanism could apply to many genetic 
traits, including, possibly, those responsible for reduction and oxidation of metals. The main 
environmental requirement that is assumed in this model for the interactions is homogeneous 
mixing of the organisms. Any environment that ensured homogeneous mixing, would be 
suitable for this model. 
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Figure 1: Schematic representation of the model assumptions. The dotted lines indicate 
bacteria entering the river from the shore. The dashed lines indicate interactions between 
bacteria which cause a transfer of the antibiotic resistant gene. The solid lines indicate 
the transfer or loss of the antibiotic resistant gene, reproduction, and death. The land 
bacteria are not adapted to the river and we assume they never reproduce. The river bacteria 
reproduce or die due to the influence of the river carrying capacity. Details are provided 
throughout the main text. 
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Parameter Description 



Dimension 



t time from head of the river 

Rs non-AR river bacteria 

Rj AR river bacteria 

Lg non-AR land bacteria 

Lj AR land bacteria 

TCB R s + Rj + L s + Lj 

a transmission rate of AR gene 

(3 loss rate of AR gene 

K carrying capacity of river 

r birth-death rate due to K 

Fs entry of non-AR bacteria from land 

Fj entry of AR bacteria from land 

7 death rate of land bacteria 



time 

concentration 

concentration 

concentration 

concentration 

concentration 

time -1 x concentration -1 

concentration x time -1 

concentration 

time -1 

concentration x time -1 

concentration x time -1 

time -1 



Table 1: Parameters, descriptions, and dimensions 



Parameter 



Value Reference 



a 


0.00006 


P 


100 


K 


1,000,000 


r 


0.01 


7 


0.02 



Andrup and Anderson (1999) 



Hendricks 



LimnoTech 



(1972) 



(2001) 



Table 2: Parameter values and references. Parameters without a reference were fixed from 
a reasonable range of possibilities for the simulations. 
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x 1 o 5 Non-AR River Bacteria 
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AR River Bacteria 



2.5 



40 60 80 100 

time (t) 

Non-AR Land Bacteria 
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Figure 2: Simulation of the populations of bacteria in a river with influx of bacteria from 
the land, over the time interval to 100. 
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Figure 3: Simulation of the populations of bacteria in a river with influx of bacteria from 
the land, over the time interval to 1000. 
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Figure 4: Parameter sensitivities - effect on non-antibiotic land bacteria, Ls- Solid curve 
corresponds to the default parameter values; dashed curve corresponds to an increase in 
the parameter value by a factor of 20%, while keeping all other parameters at their default 
value; dotted curve corresponds to a decrease in the parameter value by a factor of 20%, 
while keeping all other parameters at their default value. 
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Figure 5: Parameter sensitivities - effect on non-antibiotic river bacteria, Rs- Solid curve 
corresponds to the default parameter values; dashed curve corresponds to an increase in 
the parameter value by a factor of 20%, while keeping all other parameters at their default 
value; dotted curve corresponds to a decrease in the parameter value by a factor of 20%, 
while keeping all other parameters at their default value. 



18 



